Transient Shear Banding in a Simple Yield Stress Fluid 



Thibaut Divoux, 1 David Tamarii, 1 Catherine Barentin, 2 and Sebastien Manneville 1 

1 Universite de Lyon, Laboratoire de Physique, Ecole Normale Superieure de Lyon, 
CNRS UMR 5672, 46 Allee d'ltalie, 69364 Lyon cedex 07, France. 
2 Laboratoire de Physique de la Matiere Condensee et Nana structures, Universite de Lyon, 
CNRS UMR 5586, 43 Boulevard du 11 Novembre 1918, 69622 Villeurbanne cedex, France. 

(Dated: March 1, 2010) 

We report a large set of experimental data which demonstrates that a simple yield stress fluid, 
i.e. which does not present aging or thixotropy, exhibits transient shear banding before reaching a 
steady state characterized by a homogeneous, linear velocity profile. The duration of the transient 
regime decreases as a power law with the applied shear rate y. This power law behavior, observed 
here in carbopol dispersions, does not depend on the gap width and on the boundary conditions for 
a given sample preparation. For 7 < 0.1 s _1 , heterogeneous flows could be observed for as long as 
10 s s. These local dynamics account for the ultraslow stress relaxation observed at low shear rates. 

PACS numbers: 83.60.La, 83.50.Ax, 83.50.Rp 



Yield stress fluids (YSF) encompass a large amount 
of everyday-life and industrial complex fluids ranging 
from hair gels, cosmetic creams, or toothpastes to many 
food products or even concrete and drilling muds. They 
respond elastically below a certain stress threshold er c , 
known as the yield stress, whereas they flow as liquids if 
stressed above a c . The behavior of YSF still raises many 
fundamental and practical issues, in particular in the 
vincinity of this solid-to-fluid transition, often referred 
to as the unjamming transition [l[. Recently it was pro- 
posed to categorize YSF into two groups: (i) those with 
aging and memory effects that present thixotropic proper- 
ties (clay suspensions, colloidal gels, adhesive emulsions) 
and (ii) those for which such effects can be neglected, 
which are called simple YSF (foams, carbopol "gels", 
nonadhesive emulsions) @, Q ■ Since they present glassy- 
like features, thixotropic materials have received a lot 
more of attention in the past decade [iHU . Their yielding 
properties result from a competition between aging and 
shear-induced rejuvenation which leads to a shear band- 
ing phenomenon involving a critical shear rate j c 
Indeed, for an imposed shear rate 7 > j c , the whole fluid 
flows, while for 7 < j c , one observes the coexistence of 
a solid region with a vanishing local shear rate 7i oc = 
and a fluidized region sheared at 7i oc = j c . Besides, var- 
ious models in terms of hidden local degrees of freedom 
such as the fluidity [ll| or the structural "A parameter" 
d, 0] account for the corresponding steady-state velocity 
profiles. On the other hand, simple YSF show a contin- 
uous yielding transition in the sense that the local shear 
rate can take any nonzero value. In that case, the flow 
may appear as localized only due to stress gradients that 
result from the strong curvature of wide-gap cylindrical 
cells: in small-gap geometries, steady_-state flow profiles 
are homogeneous above yielding [lol . 12- T^j . Nonetheless, 
lately, subtle effects have been reported that somewhat 
blur the above distinction. First, non-thixotropic nonad- 
hesive emulsions, for which magnetic resonance velocime- 



try coupled to wide-gap rheology shows simple yielding 
behavior 12[ , were found to present shear banding when 
confined into microchannels 16|. Second, thixotropic 



Laponite suspensions, that exhibit shear banding in a 
rough geometry, were shown to slowly evolve towards 
homogeneous flow in the presence of smooth boundaries 
17]. These findings evidence the influence of confine- 



ment and boundary conditions on the yielding scenario 
and call for a systematic investigation of the spatiotem- 
poral dynamics of YSF. 

In this Letter we revisit the case of a simple YSF, 
namely a carbopol "gel" . Indeed, in spite of a huge body 
of rheological work, no thorough temporally and spatially 
resolved study has ever been conducted on simple YSF. 
In particular, it is generally assumed that the flow re- 
mains homogeneous during its evolution towards steady 
state. Here shear-rate controlled experiments show that 
(i) the transient regime from solidlike to liquidlike be- 
havior involves a shear-banded flow; (ii) the duration r/ 
of this transient regime decreases as a power law of the 
imposed shear rate: 77 cx 7~"; and (iii) the flow always 
turns out to be homogeneous after this transient shear 
banding. This fluidization scenario is very robust and 
the exponent a does not depend on the gap width or 
on the boundary conditions. Our results not only con- 
firm the idea that one cannot observe steady-state shear 
banding without thixotropy but also demonstrate that 
a simple YSF can present shear banding in a transient 
regime that can last surprisingly as long as 10 5 s. 

Experimental set-up and protocol - Our working sys- 
tem is made of carbopol ETD 2050 dispersed in wa- 
ter at 1% wt. This system is referred to as a carbopol 
"gel" in the literature although its microstructure shows 
a honeycomb-like assembly of swollen polymer regions 
with typical size 5-20 ^m [3, Carbopol "gels" 

are known to be non-aging, non-thixotropic simple YSF 
[l3l I20I ] . Here our samples are seeded with micronsized 
glass spheres 2l| in order to use ultrasonic speckle ve- 
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FIG. 1: The flow curve, shear stress a vs the shear rate 7, 
obtained by decreasing 7 from 100 to 0.01 s _1 (o) and then 
increasing 7 from 0.01 to 100 s _1 (•) with a waiting time of 
100 s per point. Note the slight hysteresis for 7 < 1 s _1 which 
is the signature of a significant slowing down of the system 
dynamics near the yielding transition (see text). The gray line 
is the best fit of the decreasing part by the Herschel-Bulkley 
model: a — a c + Aj n , with a c = 26.9 Pa, n = 0.55 and 
A = 9.8 Pa.s _n . Inset: zoom over the low shear rate regime. 
Both walls have a roughness of 60 [im. 



locimetry (USV) [22| simultaneously to standard Theolog- 
ical measurements with an Anton Paar MCR301 rheome- 
ter. The addition of such acoustic contrast agents only 
slightly stiffens the system (by at most 25%) and The- 
ological measurements performed without glass spheres 
show the exact same trends as in Figs. Q] and [2] 

The carbopol gel is poured into the gap of a rough 
Couette cell (rotating inner cylinder radius 23.9 mm, 
gap width e = 1.1 mm and height 28 mm). The sur- 
face roughness of 60 /xm obtained by gluing sandpaper 
to both walls was chosen to be of the order of the mi- 
crostructure size in order to avoid wall slip. To ensure 
that the strain accumulated during loading into the cell 
has no influence, the sample is systematically presheared 
for 1 min at +1000 s _1 and for 1 min at -1000 s — 1 . We 
then check that a reproducible initial state is reached by 
measuring the viscoelastic moduli G = 100 ± 10 Pa and 
G" = 12 ± 2 Pa for 2 min. Finally the sample is left at 
rest for 1 min before starting the experiment. 

Constitutive behavior - Figure [T] shows that the 
Herschel-Bulkley model accounts perfectly for the flow 
curve of our carbopol system in good agreement with 
3]. In particular, the absence of 
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previous works 13|. 
a stress plateau in the flow curve at low shear rates in a 
rough geometry is usually taken as the signature of sim- 
ple yielding behavior 0, HU . However we observe a slight 
hysteresis between decreasing and increasing shear-rate 
sweeps for 7 < 1 s" 1 (see inset of Fig. [IJ. Such an effect 
is also noticeable but not discussed in previous reports 
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FIG. 2: The shear-stress a vs time t for a shear rate 7 = 
0.7 s _1 applied at t = 0. Insets: velocity profiles v(r), where 
r is the distance to the rotor, at different times [symbol, time 
(s)]: (□, 27); (A, 1013); (©, 1730); (y 2250); (., 6778); 
(o, 9413). The arrow indicates the rotor velocity vo- Both 
walls have a roughness of 60 /im. The two-step relaxation of 
a(t) is strongly correlated to the evolution of v(r) (see text). 



I3I 2fjj . We checked that the size of the hysteresis cycle 
decreases as the waiting time is increased from 100 s to 
500 s per point. We shall see below that it may take a 
surprisingly long time to reach steady state at low shear 
rates. Therefore we will attribute the slight hysteresis in 
the flow curve to a critical slowing down of the system 
dynamics close to yielding and not to any memory effect. 

Start-up experiments - The above sweep test shows 
that a better method to probe yielding is to impose a 
constant shear rate 7 over a long duration (> 500 s) and 
to record the shear stress o~(t). Meanwhile we measure 
velocity profiles at about 15 mm from the cell bottom 
through USV [22j. Such a combined technique gives ac- 
cess to both global and local dynamics of the fluidization 
process. As shown in Fig. [2] for 7 = 0.7 s _1 , the shear 
stress follows a two-step relaxation which can be ana- 
lyzed in light of the velocity profiles: (i) for t < 2000 s, a 
shear band is present at the rotor side and slowly expands 
across the gap [Fig. |2fa)-(c)] and (ii) for t > 2000 s, 
the velocity profiles becomes linear within a few 100 s 
[Fig. HJd)] and the flow remains homogeneous and sta- 
tionary up to experimental uncertainty for the rest of 
the experiment [Fig. [2£e)]. The "fast" fluidization at 
t ~ 2000 s is better seen in Fig. Eta) where the width 
5 of the shear band is plotted against time. Comparing 
Fig. [2] and Fig. Eta), we interpret the plummet of a{t) 
around t = 2000 s as the rheological signature of the end 
of shear-banding. 

The velocity profiles of Fig. [5Ja),(b) clearly point to 
the presence of strong wall slip at the rotor even though 
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FIG. 3: Analysis of the velocity profiles recorded every 18 s 
during the experiment at j — 0.7 s _1 shown in Fig. [2] 
(a) Width 5 of the fluidized shear band normalized by the 
gap width e vs time. The shear band invades the whole gap 
over a characteristic time Tf. (b) Evolution of the effective 
shear rate 7 e ff normalized by the imposed shear rate 7 vs time, 
(c) Local shear rates 7i oc close to the walls vs time, (d) Slip 
velocities v s normalized by the rotor velocity no vs time. In 
(c) and (d), gray (black resp.) symbols correspond to the ro- 
tor (stator resp.). The errorbars indicate the uncertainty due 
to our spatial resolution of 40 /im and to the procedure used 
to fit the velocity profiles [23 ] . 



rough surfaces were used. The complete data analysis 
23J presented in Fig. |3][b)-(d) shows that smaller but 
noticeable slip is also present at the stator and that slip 
velocities as well as the local shear rate close to the rotor 
decrease concomittantly to the growth of the fluidized 
band [Fig. (3fc),(d)]. Wall slip becomes completely negli- 
gible only when homogeneous flow is reached [Fig. EJb)]. 

The above experiments and analysis were repeated for 
shear rates ranging from 0.01 to 20 s _1 and for vari- 
ous gap widths e = 1.1—3 mm and wall roughnesses 
("smooth" raw Plexiglas with a roughness of about 15 nm 
or 60 /im rough sandpaper). The scenario described in 
Figs. [2] and [3] was found to be very robust. Note however 
that, depending on the applied shear rate, the bump in 
cr(t) does not always show as clearly as in Fig. [2j for 



7 > 2 s , it gets buried into the initial decay while for 
7 < 0.2 s _1 the stress jump becomes very small. There- 
fore local velocity measurements are essential to assess 
the nature of the stress relaxation towards steady state. 

In all cases a shear band is observed transiently over 
a "fluidization time" Tf defined as the time when homo- 
geneous flow is reached as indicated in Fig. \3l&). The 
width of the fluidized zone keeps a similar evolution as 7 
is varied, although the initial value of S seems to depend 
on 7 [Fig. BJa)]. As a major result of this Letter, we 
observe that 77 decreases as a power law for increasing 
shear rates [Fig. HJb)]. For 7 < 1 s _1 , r/ is larger than 
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FIG. 4: (a) Width 5 of the shear band normalized by 
the gap width e vs time for different applied shear rates 
7 = 3.0, 2.0, 1.5, 1.0, 0.7 (gray), 0.2, 0.08 s" 1 from left to right. 
Both walls have a roughness of 60 /im. The time interval be- 
tween two successive velocity profiles was set to be at least 
1 s and inversely proportional to 7 in order to focus on the 
long-time relaxation. For each curve, the fluidization time 17 
is extracted as indicated in Fig. [3ja). (b) Fluidization time 
Tf vs the applied shear rate 7 for various gap widths and 
boundary conditions (BC). The gray line is a power-law fit 
Tf = /3/j a with a = 2.3 ± 0.1 and f3 = 272 ± 11 s 1 "". 



100 s and fluidization times as long as 10 5 s were mea- 
sured for 7 < 0.1 s _1 . Such ultraslow dynamics at low 
shear rates undoubtedly account for the small hysteresis 
in the flow curve of Fig. [T] as well as for the stress re- 
laxations reported in Fig. [5] and in 13j. The absence of 
any clear dependence on the gap width as well as the fact 
that flow curves and stress relaxations very close to those 
of Figs. [1] and [2] were measured in both cone-and-plate 
and plate-plate geometries (data not shown) imply that 
the fluidization process does not significantly depend on 
geometry-induced stress and strain heterogeneities. 

Discussion - As already implied by recent works 
T3I 2oT | , we confirm through local velocity measurements 
in small-gap Couette cells that steady-state velocity pro- 
files are homogeneous for a simple YSF. More impor- 
tantly our study brings about a number of findings that 
contradict the simple picture of a fast, continuous transi- 
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tion from solidlike to liquidlike behavior under a constant 
shear rate. A central result is the scaling law t/ oc 7~ q 
for the fluidization time. From the experimental point of 
view, it means that it is crucial to check whether the flow 
has indeed reached a steady-state, especially at low shear 
rates, before drawing conclusions from standard rheolog- 
ical data such as the flow curve or the stress response. 
Note that the stress in Fig. [2] has not fully relaxed after 
10 4 s even though the velocity profiles measured roughly 
at middle height of the Couette cell are perfectly station- 
ary. This suggests that fluidization may not occur in a 
fully homogeneous way along the vorticity direction due 
to end effects at the cell bottom or at the free surface. 

The above fluidization is reminiscent of the slow frag- 
mentation and erosion process triggered by slippery walls 
in thixotropic Laponite suspensions [l7j |. for which the 
characteristic time was also shown to scale as a power 
law with a similar exponent but with a critical shear 
rate 7* below which fluidization is never observed: rt oc 
(7 — 7*) -2 ' 1 [HI- Together with the recent simulation 
results on sheared attractive particles showing transient 
shear banding 27] and earlier observations of ultraslow 
transients in concentrated emulsions [28j . the present 
findings suggest that a generic yielding mechanism in- 
volving fluidization times with a critical-like behavior 
may be at play and go beyond the distinction between 
simple and thixotropic YSF based on the steady-state. 

Providing a model for such a mechanism, in partic- 
ular for the exponent a, appears as a major challenge 
for future work. Although a is very robust for a given 
batch as shown in Fig. @J it should be mentioned that 
values of a ~ 3 were obtained for a different sample 
preparation, hinting that a is a material parameter that 
may vary from one simple YSF to the other. Up to 
now theoretical approaches have focused on steady states 
@j E j 24 [ or on the viscosity bifurcation under imposed 
shear stress [j| [25[ so that a comparison of our experi- 
ments with existing models is not yet available. Still the 
most striking feature of the transient regime is the pres- 
ence of shear banded flows with strong wall slip in spite 
of rough boundary conditions that are usually assumed 
to prevent slip phenomena. This observation suggests a 
yielding mechanism where the initially solidlike material 
suddenly fails at the inner wall independently of bound- 
ary conditions and leaves a strongly sheared lubrication 
layer at the rotor which then fluidizes the solidlike sam- 
ple. Moreover the sudden acceleration of S(t) around 
Tf [Fig. |U[a)] an d the absence of gap dependence rule 
out propagative or diffusive processes but rather points 
to an avalanche-like behavior although no thixotropy is 
present. Such ingredients should be included in future 
tentative models for the observed dynamics. 

We thank the Anton Paar company for lending us an 
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bopol and E. Bertin, L. Bocquet, A. Colin, P. Chaudhuri 
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